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Abstract 

^f^ • In this work we solve thermo-hydrodynamical equations considering a two dimensional electron system in the integer 

quantum Hall regime, to calculate the spatial distribution of the local electron temperature. We start from the self- 
consistently calculated electrostatic and electrochemical potentials in equilibrium. Next, by imposing an external 
f— I . current, we investigate the variations of the electron temperature in the linear-response regime. Here a local relation 

C/3 ■ between the electron density and conductivity tensor elements is assumed. Following the Ohm's law we obtain local 

current densities and by implementing the results of the thermo-hydrodynamical theory, calculate the local electron 
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temperature. We observe that the local electron temperature strongly depends on the formation of compressible and 
incompressible strips. 
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1. Introduction 



> 

"^ ' Quantum Hall effects (QHE's) are remarkable phenomena observed at two-dimensional electron systems (2DESs), 

t^^ I where the longitudinal resistance /?„ vanishes in the low current regime while the Hall resistance R„. is quantized to 

• . one over integer or fractional multiples of h/e^, where e is the elementary charge and h is the Planck's constant |lll|2|]. 

^^ ' If the current (density) is increased up to a critical value, R^^ increases by several orders of rnagnitude within a narrow 



range of the current amplitude and this observation is called the breakdown of the QHE 



Tfl . To investigate 



the microscopic nature of the breakdown phenomena one essentially has to know the current distribution, i.e. where 
p\, . the Hall potential drops across the Hall bar This information is provided by local probe experiments, where the 

b ■ nri 

C^ ■ spatial distribution of the induced (Hall) potential is measured by scanning force microscope techniques BSl |9|]. It 

is observed that, the spatial distribution of the Hall potential strongly depends on the appUed magnetic field B, such 
that if the system is out of the plateau regime Hall potential varies linearly across the sample resembling the classical 
Hall effect. In contrast, within the plateaus, the potential presents strong non-linear distributions: At the high B field 
side of the plateau the potential drop is spread all across the sample indicating that the current is earned at the bulk, 
meanwhile at the low field side one observes sharp variations at the opposing edges, whereas the potential is constant 
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at the bulk. This behavior is attributed to edge state transport. The spatial variation of the potential drop as a function 
of B field shows that, while decreasing the field the edge states move towards the physical ends of the sample until the 
plateau disappears. A simple calculation of the spatial distribution of the edge states, following the original work of 
Chklovskii et al, shows that the position of the potential drop almost perfectly matches with the position of the edge 
states. The theoretical work takes into account electron-electron interactions within a Thomas-Fermi approximation 
and provides estimations of the widths and the positions of the compressible and the incompressible strips, for a fixed 
depletion length. There, it is assumed that the 2DES is separated into compressible (where the Fermi energy is pinned 
to one of the Landau levels) and incompressible strips (where Fermi energy falls into Landau gap). 

The spatial distribution of the current density and the local electron temperature at quantizing magnetic fields at 



ternperatu: 

itIIioIU] 



low temperatures in 2DES, are extensively investigated in recent years |17|, [10|, [Hi] . Giiven and Gerhardts calculated 
the spatial distribution of the current density jnaix,y) and position-dependent electrochemical potential jUec('<^,3') uti- 
lizing a local version of the Ohm's law, starting from the electrostatic quantities obtained within a self-consistent 
Thomas-Fermi-Poisson approximation (TFPA) to include classical electron-electron (Hartree) interaction 1I12I1 . This 
approach is known as the screening theory of the integer quantized Hall effect and will be briefly discussed above. In 
screening theory, heating effects were neglected hence the electron temperature T^ is assumed to be uniform through 
the system and equals to the lattice temperature Tl. Subsequently Akera and his co-workers implemented the above 
mentioned calculation scheme, however, also by solving the thermo-hydrodynamical equations they were able to de- 
scribe dissipative effects and thereby can obtain the spatial distribution of the electron temperature. They showed 
that, the current carrying incompressible strips strongly affect the electron temperature. The electrostatic boundary 
conditions imposed in this work are viable mainly for an homogeneous system, however, becomes questionable when 
considering realistic samples. In any case, this approach opens a new window to investigate the heating effects at a 
microscopic level and provides reasonable estimates when compared to experiments. 

In this paper, we investigate the spatial dependence of the electron temperature in quantum Hall systems (QHS), 
where compressible (CS) and incompressible strips (IS) are present as a result of interactions Ill3[ll4[ll5ll . We study the 
spatial distribution of the local electron temperature employing a theory of thermo-hydrodynamics in linear response 
regime, in which the current-induced change of the electrostatic potential equals that of the electrochemical potential, 
so the induced current does not influence the electron distribution lll6ll . The theory of thermo-hydrodynamics is bound 
by conservation of electron number and thermal flux densities, foUowing Akera and his co-workers OITJ llSll . Our 
investigation differs from these pioneering works, as we impose realistic boundary conditions and therefore describe 
experimental systems accurately. This enables us to predict an unexpected behavior that can be observed in the QHS 
based Aharonov-Bohm interference experiments, which is dictated by the local electron temperature variations. 



2. Model 

2.1. Basics of the screening theory 

We consider a 2DES in the z-0 plane, which is subjected to a perpendicular magnetic field B = (0, 0, B), together 
with a translational invariance in x direction. 

We assume that the electrons are confined by the background potential Vhgiy) generated due to ionized donors, 
which are distributed uniformly in the xy plane. The local electron (number) density is described by n^\(y). To describe 
the experimental geometries, we impose boundary conditions such that two metalhc gates reside at the physical edges, 
following Chklovskii et al. The eff'ective potential within the semi-classical approximation can be written as 

y(y) = VbgW + VhW (1) 

with the confinement potential 

„ / /v\2 „ 7.7rc 

iQd (2) 
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where E'^ is the pinch-off energy which defines the minimum of the bare confinement potential lll9ll . Hartree potential 
follows as 

VhW = — r dyK(y,y')n,x(y'). (3) 

K J-d 

Here k is the dielectric constant (12.4 for GaAs), no is the donor density and 2d is the sample width. Kernel K{y,y ) 
solves the Poisson's equation considering the above mentioned boundary conditions 

V(<i2-/)(c/2-/2) + (P-- y'y 



Kiy,y) = In 



(y-y')d 



(4) 



Note the fact that the boundary conditions we imposed results in different kernels compared to Ref. 112011 1811 and affect 
the strength of interactions considerably. 

The electron density is, in turn, determined by the effective potential V(y) and is calculated within the Thomas- 
Fermi approximation 

n,i(y) = r dED{E)f{E + V{y) - ^i,,), (5) 

where D(E) is the density of states (DOS), f(E) is the Fermi function and //ec is the electrochemical potential lll2lll3l 



m. 

In this work, we apply the local equilibrium approximation, used commonly to describe similar systems B22I1 . In 
local equilibrium, the energy distribution of an electron is defined by the Fermi function 

/(e,/iec,re) = , (6) 

{exp[(e-yUec) Afire] + 1) 

where s is the energy and T^. is the electron temperature. In local equilibrium approximation, the lattice temperature 
Tl remains unchanged in the presence of an applied current. If an external current is imposed the electrochemical 
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potential ^(ec(r) depends on position and its gradient E = V/iec(r)/e satisfies the local Ohm's law 



p(r)j„^,(r) = E(r), 



(7) 



hence local current densities can be obtained if local resistivities are provided. We utilize the self-consistent Bom 

approximation (SCBA) to describe both the local conductivities and density of states in the presence of an external 

I I 

magnetic field, similar to Ref. [.13.1 . Note that, the dissipative current / is the integral of current densities over the 
sample, 

/= r dyi„^.Ax,y). (8) 



= I dyjn,,x(x,y). 




Figure 1: Density profile for different lattice temperatures and the magnetic field ha>c/E^ = Q.clE^p - 0.95 {E^p - 
12.68 meV) with the sample width of d - 1.6 jim and the donor density «o - 3.61 x 10"cm"^. The insets show the 
enlarged plateau regions (incompressible strips). 



While presenting our results we use the Fermi energy at the center /s^ of the 2DES as the energy scale, since this 
quantity remains constant once the sample properties (e.g. depletion length, no etc.) are fixed, whereas the cyclotron 



frequency oJc - eB/{mc) is used to denote the field strength. The magnetic length { - y/eB/h provides a length scale, 
quantifying the importance of quantum mechanical effects such as the width of the wave function. As commonly 
used, the filling factor v = Int^n^x describes the occupation of the Landau levels. If electron density exactly equals 
to the magnetic flux (number) density the Landau level is fully occupied and v is an integer This situation is called 
an incompressible state, where no states are available at the Fermi energy. Hence, within an incompressible state, 
electron density distribution is constant. It is straight forward to define the local version of the filling factor by 
v{x,y) — 2ni^n^\{x,y). 

A typical result of the above described calculation scheme is shown in FiglT] where the electron density distribution 
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is plotted considering three different lattice temperatures. We fixed the B field strength as hojc/E^p - ^dE^p - 0.95, 
where the bulk filling factor v(0) is slightly above 2. Here, Fermi energy at the bulk is calculated to be 12.68 meV. At 
the highest lattice temperature {k^Ti^jE^p - 0.06), the electronic system is in a pure compressible state since both the 
lowest and next Landau levels are partially occupied due to the fact that the thermal energy of the electrons exceed 
the Landau gap. Therefore, no incompressible strips are formed. Once lowering the temperature, we observe that 
density distribution presents plateaus (i.e. is constant within a finite spatial interval) at opposing edges, which are the 
expected incompressible strips, where Fermi energy falls within the Landau gap locally hence, screening is poor and 
electrons cannot be redistributed. For the lowest lattice temperature, the incompressible strips become even larger and 
stable. 

Once the local electron distribution is known, using the Ohm's law and SCBA to obtain local conductivities one 
can calculate the current densities within the linear response regime where the imposed current does not effect the 
electron distribution itself. In Fig. |2] we show the calculated current density distributions considering the filling factor 
distribution in Fig. [T] One can clearly see that, at the lowest lattice temperature all the imposed current is confined 
to the incompressible strips, where scattering is completely suppressed due to the lack of available states. However, 
while increasing the temperature the current starts to flow also from the bulk. From these results, one can conclude 
that the effect of temperature is decisive both determining the electron and current densities. Therefore, in the next 
step we investigate the effect of imposed current on the electron temperature via solving the thermo-hydrodynamical 
equations. 
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Figure 2: The current densities jx versus position, for different values of the lattice temperature, Tl = 0.02, 0.03 and 
0.06 Et/kB at fixed magnetic field, Q.clE% - 0.95. Sample parameters are d - 1.6 /im and iiq - 3.61 x 10"cm"^. 



2.2. Thermo-hydrodynamical equations 

We consider two hydrodynamic equations and assume that the electron number and the total energy of the system 
at hand is conserved. The conservation of the electron number is given by 






-^ = -V.J,,,, (9) 



where the number flux density is j„., . The energy conservation is formulated by 

§ = -V.J.-Pl, (10) 

ot 

with the energy flux density j^ and the energy loss per unit area f l due to the heat transfer between electrons and 
phonons |22l]. The time evolution of the entropy density s is derived by using Eqs. (|9]l, ( fTOl i. and by the fundamental 
thermodynamical equation 

Teds - de - fiecdrieh (11) 

that yields 

T'e ^ = -V.j, - V;Uec.j„, - Pl, (12) 

where the thermal flux density j^ is described by 

jf/ = je - /"ecjn,, ■ (13) 

Utilizing the above equations that describe the drift and hopping processes, the number flux densities j^^, and the 
thermal flux densities j^ can be written as 11221] . 

j„Jx) = -0,/Zec + LllVyV - Lfj;'V,Je, (14) 

j„,{y) = -L\]y,V - LllVyfiee - Lfj;'VJe, (15) 

j,{x) = -Ll?.V,^,, + Kl^VyV - Lfj;'T/,Te, (16) 

jgiy) = -<Jv,y - 4^v_,juec - Lfj;'VyTe. (17) 



|2o|. 



Here, we use the transport coefficients (i.e. Ljt and Kji;, where {j, k] denote x and y) given by Kanamaru et al \ 

As we assume translation invariance in the x direction, the electron temperature Tg and the electrochemical poten- 
tial jUec do not depend on x. Taking into account these boundary conditions, conservation equations reduce to 

A./„,(3') = (\y\<d), (18) 

V,(A;,(y)) + e£,j,„(x) + Pl = 0, (19) 



where 



Aj«d = j«d -Jnd> (20) 






(21) 



Here j^^^ and j^'' are the number flux density and the thermal flux density in equilibrium ll20ll . 

One can see that, if the spatial distributions of the electrostatic and electrochemical potentials are known the above 
equations can be solved to obtain the local electron temperatures. In order to perform calculations for a finite width 
sample, we consider a hard wall potential {V(y) - oo for \y\ > d). Therefore flux is absent outside of the sample. 
We integrate Eqs. (18) and (19) within a narrow region on the boundary to obtain the boundary conditions. These 
conditions do not influence the distribution of electrons within the sample. 

In the following we present results, where the electron temperature is obtained starting from the results of the 
screening theory. Note that, the imposed current is assumed to be sufficiently small to guarantee that system is still 
within the linear response regime, hence, both the induced potential and the variations of the local electron temperature 
do not influence the electronic distribution. 
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Figure 3: Calculated the local electron temperature deviation Are versus position, for different values of the lattice 
temperature, Tl - 0.02,0.03 and 0.06 E^/kB- Sample parameters are d - 1.6 fim, no - 3.61 x 10"cm"^ and 
Q.c/El = 0.95. 



The local electron temperature variation as a function of spatial coordinate is shown in Fig|3] where the electron 
and current density profiles are shown in Figs. [T]and|2] respectively. We observe that, assuming a higher lattice 
temperature results in small variations at the local electron temperature. This is expected, since at higher temperatures 
the incompressible strips are not well developed, hence, current is spread all over the sample. Therefore, heating 
effects take place almost at each location. Most interestingly, one observes that one side of the sample heats up, 
whereas the opposing edge is cooled down. This effect is nothing but the Peltier effect in which cooling occurs 
at the junction of two dissimilar thermoelectric materials when an electric flows through the junction. Peltier heat 
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is proportional to the current and changes sign if the current direction is reversed Il23[ |24|| . Although there is no 
net charge transfer in the y direction the thermal flux is transferred by the conservation of energy. Lowering the 
lattice temperature yields stronger incompressible strips, hence, current is more confined to these regions. As a direct 
consequence, the local electron temperature starts to vary stronger 

Fig. ID shows the local electron temperature deviation ATf. as a function of the magnetic field with constant lattice 
temperature Tl. It is seen that AT^ shows spatially antisymmetric behavior that is obviously due to thermal flux 
transfer in the y direction. Note also that, the variation of the local electron temperature strongly depends on magnetic 
field. Since, the widths of the incompressible strips decreases by decreasing the magnetic field and the strips move 
closer to the sample edge. 
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Figure 4: Calculated the local electron temperature deviation ATg versus position, for different values of the magnetic 
field, O,/^^ =0.74, 0.8 and 0.84. Sample parameters are c/ = l.6fim,no = 3.61 x lO^cm^^ and Tl = 0.02 E^/Atb. 
The inset shows the current density profile. 



3. Conclusion 



In this paper, we calculated the local electron temperature and the current density profiles considering a quan- 
tum Hall system. First the position dependent electrostatic and electrochemical potentials are calculated within the 
Thomas-Fermi-Poisson approximation, self-consistently. Next, we used these quantities to calculate the local electron 
temperature in the linear response regime, utilizing the energy and particle conservation when solving the thermo- 
hydrodynamical equations. We found that, the current carried by the incompressible strips heats the electron system 
locally. Hence at low temperatures and at the low field side of the quantized Hall plateau the heating is mainly local, 
whereas at the high field end of the plateau heating effect is spread all over the sample. This situation is altered at 
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high lattice temperatures, since there are no well developed incompressible strips. Interestingly, we observed that one 
side of the sample heats up, whereas the opposite side is cooled down. From these observations, we expect that the 
edge-state transport is highly sensitive to formation of the incompressible strips. This leads us to conclude that, even 
at low currents local heating effects may become important. As an interesting example, we claim that the uncovered 
dephasing processes at the Aharonov-Bohm type interference experiments would yield a quantitative test for heating 
effects. We expect that, at low temperatures and at the low-field side of the plateaus the visibility of the interference 
signal is strongly influenced by the local heating effects. A simple test would be to compare the fading of the visibility 
as a function of applied current, measured at sharp edge defined (e.g. by trench gating) samples and soft edge defined 
(e.g. gate) samples. Since, the incompressible strips are more stable at gate defined samples, we predict that, heating 
effects will be more effective at sharp edge samples. Therefore, the visibility should fade away at the sharp edge 
samples faster compared to gate defined samples. 
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